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DETERMINATION OF THE ELLIPTICITY OF THE EARTH’S 
EQUATOR FROM OBSERVATIONS ON THE DRIFT 
OF THE SYNCOM II SATELLITE 


t r* 9 


by 

Carl A. Wagner 
Goddard Space Flight Center 

SUMMARY 

An original, simple theory is presented for the radial and longitudinal drift 
regime, in a triaxial earth- gravity field, of an inclined 24- hour satellite with a 
near- circular orbit. The drift regime equations for an equatorial satellite (derived 
by more complex perturbation methods) have been known for at least 2 years. This 
new theory shows that the "inclined orbit" regime is the same as the equatorial, 
modified only by an "inclination factor." 

The theory is closely validated by two numerically integrated particle drifts of 
about 3 months duration each, starting with the elements of Syncom n for epochs 
26.709 August 1963 and 10.000 December 1963. The particle program included best 
estimates of sun and moon gravity, longitude- independent (zonal) earth gravity 
through fourth order, as well as triaxial earth gravity (associated with equatorial 
ellipticity). 

On the basis of this validation, the actual drift of Syncom n over Brazil, as 
derived at Goddard Space Flight Center from range and range-rate radar and Mini- 
track observations covering a 7- month period in 1963-1964, is reduced by the theory 
to yield the following two parameters of the earth's equatorial ellipticity: 

J 22 = - (1.70 + 0.05) x 10~ 6 , 

(representing a 65 ± 2 meter difference between major and minor equatorial radii); 
and 

\ 22 = - (19 ± 6)° 

(locating the geographic longitude of the major equatorial axis). These results show 
a somewhat stronger and west- shifted equatorial ellipticity than recent geodetic 
investigations in 1963-1964 indicate. 
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DETERMINATION OF THE ELLIPTICITY OF THE EARTH’S 
EQUATOR FROM OBSERVATIONS ON THE DRIFT 
OF THE SYNCOM II SATELLITE 


by 

Carl A. Wagner 
Goddard Space Flight Center 


INTRODUCTION 

Almost as soon as Syncom n, the world's first operational 24-hour satellite, was left alone to 
drift free in the gravity fields of the earth, sun, and moon high above the earth, it began to deviate 
noticeably from the predicted path. The principal observed effect, unexplained in the orbit- 
determination theory used for this satellite (which included all well-known solar, lunar, and earth- 
gravity effects), was a small but continual daily change of the semimajor axis. As a result of this 
secular drift anomaly, the elements for Syncom n (1963 31A) have undergone revision about every 
week since mid- s umm er 1963, to allow adequate tracking acquisition for the succeeding week. 

The investigation reported herein was undertaken primarily as an attempt to explain as fully 
as possible the reasons for, and all the observed features of, this secular drift. Many possible 
causes were examined for order-of- magnitude effect and were rejected; these included magnetic 
field interactions, micrometeorite collision drag, and solar wind and radiation perturbations.* 
Outgassing from the satellite during the free drift periods starting in mid- August 1963 is believed 
to have been negligible. To have contributed to the observed steady secular increase of the semi- 
major axis, residual outgassing from the spin- stabilized satellite would have had to occur reg- 
ularly at selective times in the rotation cycle— a highly unlikely series of events over many months 
of drift. 

It was concluded that longitude-dependent earth gravity was the most likely cause of the long- 
term accelerated free drift of Syncom n. Significant longitude-dependent earth gravity will cause 
detectable perturbations of satellite orbits which can be appreciable over fairly short periods of 
time in special "Resonance" cases such as the 24-hour satellite (see Reference 1, for example). 


•See Appendix F in Goddard Space Flight Center Document X-621-64-90, April 1964, by C. A. Wagner: ’Determination of the Ttiaxiality 
of the Earth from Observations on the Drift of the Syncom II Satellite.” 
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Summary of Previous Longitude-Dependent Earth-Gravity Investigations 


The question of the existence and extent of the longitude dependence of the earth's gravity 
field has concerned geodesists since the early years of this century (see Reference 2, for example). 
The existence of a longitude- dependent field implies the existence of inhomogeneities and states 
of stress within the earth which are of fundamental importance to all dynamical theories of the 
earth's interior. 

Table Al, Appendix A, summarizes 40 previous reductions of gravimetric, astrogeodetic, and 
satellite gravity data reporting longitude terms in the earth's external gravity field. It is quickly 
seen that, even though the more recent gravity reductions show closer agreement in term-by- 
term detail than the older ones, exact knowledge of tesseral earth gravity is still in its infancy— 
basically, because of the scarcity of good data over more than a small portion of the external 
field. Surface gravity data, prior to the recent use of accurate sea-going gravimeters, have been 
very scattered and often of dubious quality. Uncertainties in surface- station position with respect 
to the center of the earth's mass of the order of only 10 meters can cause serious errors in tes- 
seral gravity reductions, since this is the order of magnitude of the geoid deviations caused by 
longitude-dependent gravity. In addition, scattered surface gravity data, unless smoothed with 
great care, tend to overemphasize the higher order tesseral terms of greater spatial frequency, 
because the gravimeter will be sensitive to the gravity distortion of even close-by mass anomalies. 
Satellite gravity reductions, inherently more insensitive to the higher order terms, have suffered 
most to date from being based on too limited a number of medium- altitude medium- inclination 
satellites. The chief result of this poor field sampling has been to make it difficult to separate the 
effects of individual tesseral terms of the same spatial frequency. 

The situation does not seem as bad when the constants of the most recent tesseral gravity 
reductions are taken as a set and displayed cumulatively as deviations from a mean earth geoid 
(Appendix A, Figures A2 to A8). Still, as late as July 1963, Izsak stated (Reference 3), "It might 
be some time before one can arrive at definite conclusions regarding the longitude dependence of 
the earth's gravity field." 

The use of 24- hour near- circular orbit satellites for geodetic purposes (as revealed in this 
investigation) has certain obvious advantages over lower altitude and surface data reductions. For 
the near-stationary satellite, any small-latitude symmetric earth- gravity anomaly in longitude 
will, in time, cause significant longitude drift of the ground track configuration which can be 
evaluated to high precision after a long period of observation. The great height of the 24-hour 
satellite makes it sensitive only to the very lowest orders of gravity anomalies, thus making it 
easier to separate out the effects of individual tesseral terms in the total field (see Appendix A, 
Figure Al for example). 

On the other hand, the great altitude of the 24-hour satellite has its drawbacks. The longitude 
perturbations to be sensed are extremely small at 24-hour altitudes: of the order of 10 -7 of the 
earth's principal gravity attraction. The observation time for an accurate determination of this 
perturbation field must be long. Furthermore, the effect on the distant satellite of sun and moon 
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gravity over these long observation times cannot always be safely ignored, as it generally can for 
reductions from lower altitude satellite data. 

In summary, it may be said that all earth satellites orbiting well above the atmosphere and 
below a level where sun and moon perturbations become too severe can be potentially useful in 
helping to clarify knowledge of the earth's actual gravity field. Twenty- four hour satellites appear 
ideally suited to defining individual components of this field to at least third order. 

Purpose of Study 

The author is not aware of an adequate presentation in the literature of the drift theory for an 
inclined 24-hour orbit satellite (such as Syncom n) in a longitude- dependent earth gravity field. 

The purpose of this report, then, is threefold: 

1. To present a theory for the drift of a 24-hour inclined satellite in a longitude- dependent 
earth- gravity field; 

2. To interpret the actual drift of Syncom II (as assessed by a set of revised elements for the 
satellite) on the basis of this theory, thereby deriving measures of the longitude- dependent 
gravity field (i.e., shape of the geoid) which is assumed to be responsible for the observed 
drift; 

3. To calculate maximum on-board station- keeping propulsion requirements for future 24- 
hour satellites based on the drift acceleration experience of Syncom n. 

The measures of the earth's shape (geoid) derived in this report represent the first use of a 
24-hour satellite for this purpose. 


BASIC THEORY OF THE REDUCTION: ORBIT-AVERAGED DRIFT OF A 24-HOUR SATELLITE 

(Determination of the Longitude Drift and Orbit Expansion for a 24- Hour Satellite in a Near- 
Circular Inclined Orbit Affected by a Small but Persistent Tangential Perturbing Force) 

The dominant perturbations of a 24- hour equatorial satellite in a higher order earth- gravity 
field have been derived many times in the literature (References 1, 4, and 5).* In these references 
the perturbations were found by directly linearizing the equations of motion themselves and dis- 
playing the perturbed motion in appropriate geographic coordinates; no attempt was made to treat 
the drift of the inclined 24-hour satellite. 

This report departs from the rather involved and difficult-to-visualize procedure of lineari- 
zation of the equations of motion. Instead, it is shown how simply the dominant drift and orbit 


*Also: Barrett, C. C., ‘’The Perturbations of a Synchronous Satellite Resulting from the Gravitational Field of a Triaxial Earth,* GSFC 
Document X-623-62-160, September 1962. 

And: Private communication from R. H. Frick and T. B. Garber, 1962 (Rand Corp. Memo RM-2296). 
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expansion equations for the 24-hour satellite can be found by calculating the "perturbation of the 
two-body energy" of the geographically stationary satellite due to the small but persistent longitude- 
dependent earth- gravity force. This paper will not discuss in detail the limits of validity of the 
expressions derived, except in the case of the "inclination factor," which is discussed in Appendix 
B. To assess the accuracy with which these expressions predict the satellite's behavior, simu- 
lated trajectories with typical Syncom II orbit elements have been run on Goddard's particle pro- 
gram ITEM (Interplanetary Trajectory by an Encke Method). These trajectories (Appendix C) 
confirm the validity of the derived drift equations to an accuracy well within the "noise levels" in 
the orbital elements reported for Syncom II (Appendix D). The equations are essentially the same 
as those which Dr. Musen has derived (though not explicitly) by a more general but complex 
"energy perturbation" method (Reference 6). 

In Figure 1, F is a small earth-gravity perturbation force acting tangentially to an initial cir- 
cular 24-hour satellite orbit and ds is a small arc length of the satellite's path around the earth.* 

At the beginning of the dynamics, the total energy (the sum of potential and kinetic) of the satellite 
in a spherical earth- gravity field (Reference 7) is 



where [± E is the earth's gaussian gravity constant (3.986 x 10 s km 3 /sec 2 ). The energy added to the 
satellite by f per day is 


AE = ®Fds 


PROGRESSION OF 



Figure 1 — Orbit plane of a 24-hour satellite, 
looking southerly. 


- Ztt a s r , ( 2 ) 

where F = (1/2-rr)^ Fdi9. In units of force per 
unit mass, F is the orbit- averaged energy per- 
turbing force. If the orbit is purely circular, 
only a tangential perturbation force can cause a 
change in the total energy. The ITEM simu- 
lated trajectories in Appendix C and the real 
Syncom n orbits both maintain eccentricities of 
the order of 0.0001 for periods up to 100 days. 
Equation 2 assumes that the eccentricity is zero 
for the 24-hour satellite of semimajor axis a s . 

From Equation 1, the change in energy of a 
24 -hour satellite is accompanied by a change in 


•Symbols used in this report are defined in Appendix F. 
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M E Aa, 


4 ? * 

semimajor axis expressed by 


or 


AE = 


k y 


Aa . 


2 (a,) 2 AE 


**e 


(3) 


Substituting Equation 2 into 3, the change in semimajor axis of the 24-hour near-circular orbit, 
per day, is approximately given by 


Aa. 


4Wa.\ 3 F 


IN 


(4) 


From Kepler's third law, the period of a 24-hour orbit as a function of its semimajor a -rig is 


w 


V2 


(5) 


Thus, if the semimajor axis chaises by Aa s , the period change is given by 

3^(a t ) 1/2 Aa s 


AT 




V2 


( 6 ) 


Substituting Equation 4 into 6, the change in period, per day, of a 24-hour circular orbit is given by 

12-77 2 (a s ) 7/2 F 


AT 


K) 


3/2 


(7) 


The apparent net longitudinal drift rate of the 24- hour satellite’s ground track with respect to the 
surface of the earth (see Figure 3) after the first sidereal day is 


K(t - 1 sid. day) 


(AT ) 2n 


rad/sid. day.* 


( 8 ) 


The minus sign is taken in Equation 8 because a gain in period is accompanied by a decrease in 
net geographic longitude for the initially 24-hour satellite (for example, for the daily geographic 


•See Is ley, W. C., *A summary of Constants Associated with Orbital Analysis of Earth Satellites, Including the Influence of their Un- 
certainties upon Gravitational Measurements for Synchronous Satellites, * GSFC Document X-623-62-169, 1962. 
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position of the ascending equator crossing). Combining Equations 7 and 5 in 8 gives 


\(t - 1 sid. day) 


12t t 2 F 

v(>.) ! 


rad/sid. day. 


( 9 ) 


As the gain in semimajor axis is small over one day (and, in fact, small compared with a s for 
the entire libration period of the satellite in the triaxial earth field), the drift rate will continue to 
build up linearly with time initially, adding increments of Equation 9 each day. Thus, the net longi- 
tudinal drift acceleration of an initially 24-hour satellite is 


K = 


12772 f 

^e/K) 2 


rad/sid. day 2 . 


( 10 ) 


For the equatorial synchronous satellite, F is a constant at every point in the orbit equal to 
the longitude perturbation force at that equatorial position. For this satellite, the drift accelera- 
tion as a function of k (in a second-order gravity field) is simply given from Equation 10 as 


k - A 22 sin2 (\-X. 22 ) rad/sid. day 2 , 


(10a) 


where 


A 22 (equatorial) = - 72 tt 2 J 22 (R 0 y& s ) 2 rad/sid. day 2 . 


(10b) 


(See next section and Appendix A, Equation A4.) Equation 10a has been derived by a more complex 
perturbation method.* Rewriting Equation 4 as 


47r(a s ) 3 F 

a _ — length units/sid. day 


(ii) 


gives the expansion rate of the initially 24-hour near-circular satellite orbit due to a small but 
persistently acting orbit- averaged tangential perturbing force f. 


EVALUATION OF ORBIT-AVERAGED PERTURBING FORCE FOR AN 
INCLINED-ORBIT 24-HOUR SATELLITE 

While the energy- changing force on the stationary equatorial satellite is constant over a 
single orbit, this is no longer true for the inclined 24-hour satellite. Such a circular orbit satel- 
lite describes a closed, narrow "figure 8" path over the earth's surface, centered on the equator 
(see Figure 3). 


*Frick and Garber, op. cit. (See footnote, p. 3). 
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At each point on this path the tangential, energy- changing force is different, as both the longitude 
and latitude change during the daily excursion. This along-track force is now made up of contri- 
butions from both latitude and longitude gravity perturbation forces (Appendix A). It can be shown 
that the zonal gravity forces have no net daily energy effect on the circular inclined orbit. The 
net daily contributions from the latitude and longitude perturbation forces due to the earth's pre- 
sumed elliptical equator are not zero. It turns out that the net contribution from the longitude 
perturbation dominates for orbits of s mall and medium inclination. 

Figure 2 shows the position of the 24-hour satellite with respect to the earth and the celestial 
sphere. The earth-gravity perturbing forces in the radial, latitude, and longitude directions, F r , 
f^, and f^, are assumed to be acting on the satellite at s. 


NORTH 



Figure 2 — Position of a 24-hour satellite with a near-circular orbit with respect 
to the earth and the celestial sphere. 


If only the earth- gravity perturbation forces arising from the ellipticity of the earth's equator are 
considered,* Appendix A gives these forces as 


(*.)■ 


[9J 22 cos 2 <£cos 2 (X A. 22 )J 


( 12 ) 


* Wagner, C. A., '’The Gravitational Potential of a Triazial Earth,” GSFC Document X-623-62-206, October 1962. 
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^ _ Me ( R o/ a s ) 2 p , 

- : — “ |^6J 22 cos 0 sin <fi cos 2 (A. - /^ 22 ) J ’ 


F A ° 


( a s ) 2 

( R o/ a s ) ' 


^6 J 2 2 cos 0 sin 2 ^ 22 )] 


(12a) 


(13) 


As long as the orbit is nearly circular, F r will have negligible contribution to f. The contribution 


to F from F 0 is 


F ^ - F^ cos a - K sin 0 cos 0 cos a cos 2 (A. - \ 22 ) ; 


(14) 


K is a constant for a single orbit: 


Me Rq 2 

K = 6J ” KT 


(14a) 


In the right spherical triangle AN, S,L note the following trigonometric relations: 


From Equation 15a, 



tan AL 


cos i 

tan 6 

(15a) 


tan 0 


cos a = 

tan 6 

(15b) 

sin 0 - 

s in i sin 6 , 

(15c) 


cos i 


sin a = 

COS 0 

(15d) 

- tan -1 

(tan 6 cos i ) . 

(16) 


Let the geographic longitude of the satellite at the ascending node an be \ 0 . If time is counted 
from this orbital position, the geographic longitude of the 24-hour satellite at s in its near- 
circular orbit is 


A - + AL - oi e t (Figure 2) ; 


or, using Equation 16, 


A. - \ 0 + tan 1 ( tan 6 cos i ) - ai e t ; 


(17) 
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(a> e is the earth's sidereal rotation rate). For the 24-hour satellite [starting the dynamics with S 
at AN (ascending node) for convenience], e - <o e t , so that Equation 17 becomes 


A - A 0 + tan *(tani9cosi) - 6 , 

approximately. The function tan -1 (tan 6 cos i)- e is "even" about 0=0 and e = V2, with a 
period of -n, and behaves like a somewhat distorted sine function (Figure 3 and Appendix B). Call 
this function AA and note that, for i < 33 degrees, A A is always less than 5 degrees. Thus, using 
Equation 18 and assuming i is sufficiently small (i < 45 degrees proves to be a sufficient restric- 
tion on the inclination), cos 2 (a - a 22 ) for the 24-hour satellite can be approximated by 

cos 2 (A - A 22 ) = cos 2(A 0 -A 22 ) - 2AA sin 2 (A„ - A 22 ) = - cos 2y 0 + 2AA sin 2r 0 . (18a) 

(Note that y 0 = 90°+a o -a 22 fromFigure3.) Similarly, sin 2 (a- a 22 ) can be approximated by 

sin 2 (A 0 - A 22 ) + 2AA cos 2 (A 0 - A 22 ] = _ s i n 2y 0 - 2 Aa cos 2y 0 = sin2(A-A 22 ) . (18b) 

In Equations 18a and 18b, y 0 is the geographic longitude of the node of the 24-hour near- 

circular satellite orbit with respect to the minor equatorial axis (Figure 3). With these expansions 
(Equations 18a and 18b), and using 15b, Equation 14 becomes 

„ sin 4 > cos 4> tan 4> , . _ , 

F ( *> - K * (~ cos 2 ^o + 2^A sin 2y Q ) . 


Using Equation 15c in the above expression, the contribution to the perturbing force F due to 
becomes 


„ sin 2 i sin 26 , 

'm - K ' 2 v cos 2y 0 


+ 2A\sin2y Q j . (19) 

Writing AX = A\ max sin2<9, Equation 19 becomes 


• r , 0 % sin 26 

F &) - K sin 1 r cos ^0) 2 


K sin 2 i A\ max (sin 2y 0 sin 2 2(9) . (20) 


Averaging F^ over 0 < 6 < 2tt, Equation 20 gives 


j.r 

277 1 

Jo 


I 2 iAA 


F C *> = % 1 F ^) d6 = 


max f \ 

~~2 (K sin 2y 0 ) . ( 2 1) 



LONGITUDE OF 
EQUATORIAL 
MINOR AXIS 


PROGRESSION OF 
SUBSATELLITE 
POINT AWAY FROM 
ASCENDING NODE FOR 
A 24 -HOUR SATELLITE 


INSTANTANEOUS 

SUBSATELLITE 

POSITION 


EAST 


T 

EQUATOR 


Figure 3— Geographic subsatellite track of a 24-hour 
satellite in a near-circular orbit. 
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The contribution to F from F x is 

F^ = F^ cos (90° -a) = F^sina = K sin a cos <p (- sin 2y Q - 2A\ cos , (22) 

from Equations 13 and 18b. Using Equation 15d in 22 and noting that A \ ± - A^ max sin 26, as before, 
gives the contribution to F from f a as 

F C0 = K cos i [- sin 27 0 + 2A\ max sin 20COS 2y 0 ] . (23) 

Averaging f (A) over 0< 6 < 2-n, Equation 23 gives 

F (A) = K cos i sin 2y 0 . (24) 

Thus, combining the contributions of the latitude and longitude perturbations to the average 
perturbation force over a single 24-hour orbit, Equations 21 and 24 sum to produce 


F 


total 


F «>) + F <A) 


- K sin 2y 0 cos i + 


sin 2 i Ax 


(25) 


COMPLETION OF DERIVATION OF DRIFT EQUATIONS FOR A 24-HOUR SATELLITE 

Appendix B shows that 

A\ max = tan -1 (sec i) - 45 degrees . 

It is also shown there that, to a high degree of accuracy for i < 50 degrees, 


cos i + 


sin 2 

max 

5 


1 + cos 2 i 
2 


Numerically integrated orbits have shown that the drift theory for a 24-hour satellite stemming 
from Equation 25 is in error by more than 2 percent for i > 45 degrees. With this restriction in 
orbit inclination, using the above approximation for the inclination factor, we can rewrite the 
longitude drift and orbit expansion equations (10 and 11), evaluating F by Equation 25, giving 

.. _ 12tt 2 K ( cos 2 i + lN . 

7 ° " M E /(a") 2 { 2 J sln 2y 0 rad/s id. day 2 (26 ) 


(where y 0 refers to the net acceleration of the initial ascending equator crossing longitude y 0 ) and 
• - ^( a s ) 3K 


^cosM_+i^ 


sin 2y 0 length units/sid. day 


(27) 
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Substituting Equation 14a into 26 and 27 reduces these expressions to 


and 


y Q - 72tt 2 J 22 (R 0 /fe s J 2 ^ 2 jsin2y 0 radsid. day 2 


/« , \ - f COS 2 i + 

J 22 \ R o/*s ) Ro 1 2 ) sin 2 ^o len 6 th 


= - 2477 T.o R« I o I sin 2y n leneth units /sid. day 

Define a nondimens ional change of semimajor axis from a s during the drift as 


Aa 


a i - — " F - so that a 1 - 3- 


With Equation 29a, Equation 29 becomes 


aj 24 t 7 J 22 (R 0/ /& s ) ‘ 


cos 2 i + 1 


sin2y 0 , l 7 sid. day 


( 28 ) 


(29) 


(29a) 


(30) 


Define 


A 22 72m 2 J 22 (R„/a s ) 


cos 2 i + 1 


rad/sid. day 2 


With Equation 30a, Equations 28 and 30 become 


(30a) 


y 0 + A J2 sin 2ry 0 ~ 0 rad/sid. day 2 (31) 

and 

A 22 sin2y 0 

3j - 5^ = 0 , 1/sid. day (32) 

The long-term drift regime for the 24-hour inclined- orbit satellite (as described by Equations 
31, 32, and 30a) is thus the same as for the equatorial satellite* modified only by the inclination 
factor (cos 2 i + 0/2- 

Note that m in Equation 32 has dimensions of radians per sidereal day. It must be understood 
that Equation 31 describes the net daily geographic acceleration of the initially 24-hour satellite 
with respect to the earth's minor equatorial axis. Stated another way, Equation 31 describes the 
geographic drift of the entire originally stationary, figure- 8 ground track (Figure 3). Similarly, 


•See: Frick and Garber, op. cit and Barrett, op. cit. (footnote, p, 3)* 


11 



Equation 32 describes the net daily orbit- expansion rate of the 24-hour satellite. In particular, it 
is convenient to treat the motion of the ascending node of the orbit in geographic longitude as a 
reference for the entire configuration. In what follows, therefore, y will refer always to the geo- 
graphic longitude of the ascending equator crossing east of the equatorial minor axis and y 0 will 
refer to the initial geographic longitude of the ascending equator crossing east of the minor axis, at 
the start of the dynamics under consideration. Equations 31 and 32 can thus be rewritten in terms 
of the general ascending equator crossing longitude position y, to give the relevant partially un- 
coupled long-term drift and orbit- expansion differential equations for the near-24-hour near- 
circular orbit satellite: 


and 


y + A 22 sin2y - 0 rad/sid. day 2 

A 22 sin 2 y 

” 3^ - 0 , 1/sid. day 


GENERAL CONSIDERATIONS OF SOLUTIONS OF DRIFT EQUATIONS: LIBRATORY 
DRIFT REGIME OF AN INCLINED-ORBIT 24-HOUR SATELLITE 


Equation 33 can be integrated directly for the geographic drift rate by noting that 

•• _ d2 _ d y 2 _ d y 2 

j ~ dt ~ 2y dt 2d y 

Thus, Equation 33 can be separated to 


d y 2 - - 2 A 22 sin 2y dry . 

Since the variables y 2 and y are separated in Equation 35, this equation integrates to 

y 2 = A 22 cos 2y + Cj . 

With the initial condition that y = y 0 at y = y 0 , Equation 36 becomes 


(33) 

(34) 


(35) 


(36) 


y - [(>o) 2 + A 22 ( cos 2 ?~ cos 2 >o)] 1/2 ’ (37) 

giving the drift rate of the 24 -hour satellite as a function of the initial drift rate y 0 , the earth- 
gravity constant A 22 , the initial longitude east of the minor axis y 0 , and the instantaneous longi- 
tude y. Returning to the semicoupled system of Equations 33 and 34, the explicit dependence of 
the equations on the location from the minor axis and the magnitude of the equatorial ellipticity 
may be eliminated by multiplying Equation 33 by l/3w and adding the resulting equation to 34, giving 

y + 3tt = 0 . (38) 
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Equation 38 can be rewritten as 



Separation of the variables y 


377 a x - 0 - dy + 377 a x dt - dy + 377 da x . (39) 

and aj is thus achieved in Equation 39, which integrates directly to 


^ a i + y - c 2 ■ (40) 

With the initial conditions, a x = 0 when y = 0 (or the satellite is in the momentarily station- 
ary ground-track configuration), the integration constant C 2 is evaluated as 0. If y 0 in Equation 37 
is also the longitude of this initially stationary orbit, (y 0 ) 2 = 0 there, and Equation 37 in 40 yields 
for a j , the semimajor axis change from "synchronism" in the drift motion, 

(A 22 )‘ /2 ( c ° s 2r- c °s2y 0 ) 1/2 . 

ai - + 3^ • 

Equation 41 shows explicitly that the semimajor axis is bounded in long-term drift from a station- 
ary orbit. From Equation 33, since A 22 > 0, if 0 < y 0 < 90 degrees, drift proceeds toward the near- 
est longitude of the earth's equatorial minor axis (in the -y direction). If -90 degrees < y g < 0 degrees, 
then 33 shows that drift again proceeds toward the nearest minor axis longitude (in the + y direction). 
Thus, in all cases of drift from a stationary geographic configuration, cos 2 y - cos 2 y 0 is a positive 
function which has a maximum when y = 0 (when the satellite has drifted over the longitude of the 
minor axis). Thus Equation 41 gives (for the librations of a 24-hour satellite) 


a 


1 ,max 


(A„) 1/2 (1-cos 2y 0 ) 1/a 
3m 


(42) 


Again, it is noted that v in Equation 42 has units of radians per sidereal day. An absolute maxi- 
mum semimajor axis change in the drift occurs when the "synchronous" condition is established 
near the longitude of the major equatorial axis. Here, y 0 = -90 degrees, cos 2r 0 = -1, and 

a 2 (absolute maximum for a librating 24-hour satellite) = ^ • (43) 

For the constant J 22 = -1.7 x 10*® (derived in this study from long-term observations on the 
drift of the Syncom n satellite) and using the additional constants from this study (i = 33 degrees, 
a s = 42166 km, R 0 = 6378.388 km), Equation 30a gives 

A 22 = 23.2 x 10"® rad/sid. day 2 . 


Equation 43 then gives 


a 2 (absolute max. ) - 0.72 x 10 3 
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from which, by Equation 29a, 


A a (absolute max. from a “synchronous ” condition near 

the equatorial major axis, for a satellite of i = 33 degrees) = 30.7 km. 


Thus the assumption made in Equations 10 and 11, to approximate the slightly varying semimajor 
axis by a s (a constant) throughout the drift motion, appears amply justified. 

Figure 4 is a graph of Equation 41 for a x vs. y (the longitude with respect to the nearest minor 
axis location) as a function of y 0 , the longitude in the initially stationary configuration. 

Note that Equation 41 allows equal ± values for a x for each y. Suppose the satellite is initially 
at +y Q (position 1 in Figure 4) from the nearest location of the minor axis. From Equation 33, 
sin 2y 0 being positive, the satellite begins to drift west (attaining a negative drift rate) toward the 
minor axis. But, from Equation 40, since C 2 = 0, a x = - 7/377 >0; the drift therefore proceeds 
counterclockwise in Figure 4, around the central point of the minor axis and a x = 0, along the up- 
per portion of the two- valued arc determined from Equation 41. 
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Figure 4— Libra ti on with longitude of the semimajor axis of a 24-hour satellite as a 
function of the longitude of the initially stationary configuration. 






I • 

The same situation holds for the motion beginning or stemming away from the "synchronous" 
longitude at - y 0 , position 2 in Figure 4. Here sin 2y 0 is negative, and the drift proceeds at a posi- 
tive rate to the east. Again from Equation 40, as soon as the satellite leaves position 2, a 1 = -y/5n< 0, 
and the circulation continues in a counterclockwise direction. Every trajectory in the phase plane 
7 may be conveniently defined by the constant q of the "energy integral" of the drift motion 
(Equation 36). Since Equation 33 is the equation of motion defining the large-angle oscillations of 
a mathematical pendulum (in the case of the 24-hour-orbit satellite, the point of symmetry is the 
minor axis where 2y = 0), it can be expected that the general solutions in that theory (Appendix E) 
apply to a good approximation to the long-term librations of the near-synchronous satellite, pro- 
vided the drift rates are not large.* For example, in Equation 36, with a momentarily "synchronous" 
condition at y 0 being given by y 0 = 0 , the "energy constant" is evaluated as 

C 1 = - A 22 005 ^0 ' 

With this evaluation, Equation 36 becomes 

7 2 ~ A 22 COS 2 > = ' A 22 COS 2 7o * (44) 

Solving for the initially "synchronous" longitude as a function of any longitude in the drift and the 
corresponding longitude rate, Equation 44 gives 



Since y 2 /A 22 > 0 , the argument of cos 1 in Equation 45 is always less than or equal to 1. Thus, as 
long as (cos 2y - y 2 )/^22 > “ 1> then 45 will give a real solution for the momentarily "synchronous" 
longitude with respect to the minor axis. But, if (cos 2 y - 7 2 )/b 2 2 < ' 1> there will be no real momen- 
tarily "synchronous" configuration for the near- 24-hour satellite. With this energy, the world- 
circulation regime commences, corresponding to the over- the- top, high-energy regime of the 
mathematical pendulum (Reference 1). The above inequality implies that, for the commencement 
of "world circulation" for the near-24-hour satellite, 

y 2 > A 22 ( l + cos 2y) , 


or 


y 2 - 2A 22 co * 2 y * ^ 

When 2y = 0, or the satellite is over the minor axis, Equation 46 allows the maximum possible 
drift rate for a librating 24- hour satellite: 


^max, for libraticn ° «d/sid. day 


(47) 


•Author has two papers "in preparation, 1965” on this subject, to be Goddard Space Flight Center X-documents. 
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For example, using the reported value of A 22 = 23.2 X 10~ 6 rad/sid. day 2 for the inclination of 
the Syncom II satellite, Equation 47 gives 

y o " (46. 4x 10“ 6 ) 1/2 s 0. 39 degree/day . /, ft \ 

max, for libration with J ^ " - l-7xl® , i - 33 \ * ® / 

APPROXIMATIONS TO EXACT DRIFT SOLUTIONS FOR PERIODS 
VERY CLOSE TO SYNCHRONOUS 

Although a first integral to Equation 33 can be found easily (Appendix E), the second integral 
can give the drift only as a function of time in a closed form, in terms of an elliptic function. How- 
ever, for reasonably small excursions from a momentarily synchronous configuration, Appendix E 
shows that the elliptic function can be approximated to good accuracy by the first few terms of a 
Taylor series expansion in the time from synchronism At . 

Expanding the drift from the "synchronous" longitude (y - y 0 in this section) in a Taylor 
series, with respect to increments of time At from the momentarily stationary condition, 


y(At) 


.. (At) 2 

y 0 + y 0 At + y 0 — ~ 



y M 

x 0 


(At) 4 

24 


+ 


To [5] 


(At) 5 

120 


?o 720 


(At) 7 
5040 + 


[ 8 ] 


(At) 8 

40320 


(49) 


Differentiating Equation 33 six times with respect to time, it is clear that all derivatives in Equa- 
tion 49 can be written as functions of ^ 22 , y 0 , and y 0 . Noting that y(At) -y 0 = Ax. (the geographic 
longitude with respect to the ’’synchronous" configuration) and y 0 = 0, then 49 can be shown to 
reduce to the expansion 


(At) 2 


(At) 4 , 

24 + 

[( a 22 ) 3 s in 2y 0 (4 s in 2 2y 0 - 1 ) 

I (At) 6 
1 180 


2 2 >o-l)' 

1 (At) 8 

J TUSSO 


(50) 


It is apparent that, as At - 0, the higher order terms of Equation 50 become increasingly more in- 
significant to the total drift, in comparison to the terms of lower order. 

In Appendix E, the exact "elliptic integral" of motion from Equation 33 is calculated from a 
synchronous longitude of 60 degrees east of the minor axis. This calculation demonstrates that 
the simple term- inclusion-time criterion below gives an adequately converging series to the 
"exact" drift. In the actual reduction, all higher order terms in Equation 50 which are less in 
magnitude than the root mean square (rms) error of the observed Syncom II longitudes are ignored. 
The section on reduction of orbits (p. 22) in this report shows that this rms error of longitude de- 
termination for the ascending equator crossings of Syncom II from August 1963 to March 1964 has 
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been of the order of ±0.025 degree. Thus, 0.025 degree is used below in for min g the minimum- 
time - te r m- inclusion criterion for each term of Equation 50. A 22 is assumed to be 23.2 x 10 -6 
rad/sid. day 2 . 

1. For inclusion of the (At) 4 term, 


|sin4y 0 | is maximum when j 0 - ±22.5 and +67.5 degrees. 


Therefore, 


| (from the fourth-order term)| = (A 22 ) 

Solving Equation 51a for At, when |A^ max (fourth order) | = 0.025 degree, 


(51a) 


1/4 


At (min. fourth-order term inclusion) - ^0. 025 x 24/57. 3x (23.2 x 10 6 ) J 

- 66.5 sid. days from “synchronism 


2. For inclusion of the (At ) 6 term, 


sin 2j 0 (4 sin 2 2y 0 - l) 


is maximum when y 0 - ±45 degrees. 


Therefore, 


A\ (from the sixth-order term) 




(51b) 


Solving Equation 51b for At , when | A^ max (sixth order)| = 0.025 degree, 

At (min. for sixth-order term inclusion) = |^0. 025 x 60/57 . 3 x ( 23. 2 x 10“ 6 ) j 


1/6 


113 sid. days from “synchronism. 


3. For inclusion of the (At) 8 term. 


sin 4 y 0 (34 sin 2 2y 0 - l) 


is maximum when y Q - ±59 1/4 degrees. 


Therefore, 


A/ (from the eighth-order term) 


(At) 8 

212 ( a 2 2 ) 10050 


(51c) 
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Solving Equation 51c for At, when |A\ max (eighth order)| = 0.025 degree, 

At (min. for eighth-order term inclusion) = |^10080 x 0 . 025/21 . 2 x 57 . 3 x ( 23 . 2 x 10 6 ) J 

- 171 sid. days from “synchronism.” 

Similarly, expanding a 1 (t) in a Taylor series about the time of "synchronism" ( t 0 , y 0 ), 

(At ) 2 (At ) 3 

a i (to + At ) “ ( a i)o + ( a i)o (^ t ) + ( a i )o 2 + ( a i ) o 5 + • (52) 

But, from Equation 34, 


, , A 22 sin 2r0 

( a i)o - • (52a) 


Differentiating Equation 52a with respect to time gives 


2 7 0 A 22 cos 2 y 0 

(«i)o = 5T = 0 - (52b) 


since y 0 = 0. Differentiating Equation 52b with respect to time gives 


~ 4 ('Xo) 2 A 22 sin 2 yp 2y 0 A 22 cos 2y 0 - (A 22 ) 2 sin 4y 0 

( a 1 ) 0 “ 3m + 3tt ~ 3m ’ (53) 


using Equation 33. From the conventional definition of a : , (aj) 0 = 0. Equation 52 then becomes 


(A 22 sin 2y 0 ) At ( A 22 ) 2 sin 4y 0 (At ) 3 

aj (at At from "synchronism”) - ^ + ••• , (53a) 


with the results of Equations 52a, 52b, and 53 in 52. 

The section on reduction of orbits shows that the rms error of semimajor axis determination 
for Syncom II (including sun and moon "noise") is the order of ±0.5 km. Therefore, the rms error 
to be expected in is of the order of 0.5/42166 = 1.185 X 10~ 5 . Following the procedure for the 
longitude drift, 1.185 x 10' 5 is used below to determine the minimum time for the inclusion of the 
terms beyond the first on the right-hand side of Equation 53a, to insure adequate convergence of the 
infinite series for (At). 

1. For inclusion of the (At) 3 term, 


sin 4y g is maximum when y 0 - +22.5 and ±67.5 degrees. 
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Therefore, 


a, (from the third-order term of 53a) I 

I i, max y \ 

Solving Equation 54 for At, when |a 1>max | = 1.185 x 10~ 5 , 

At (min. for the third-order term inclusion in 53a) - 


From a "synchronous" configuration at 54.8 degrees west of Greenwich, on or about September 6, 
1963, Syncom n drifted to 59.2 degrees west of Greenwich on November 28, 1963, where it was 
"stopped" by the tangential firing of on-board cold-gas jets. A second free-drift period followed 
from a "synchronous" configuration at 59.2 degrees west on about November 29, 1963, to 66.3 de- 
grees west on March 18, 1964, where the on-board tangential jets were fired to speed up the west- 
ward drift. Of the 34 separate orbits calculated by the Goddard Data and Tracking Systems 
Directorate for these free-drift periods, only 7 fell outside the minimum 66- day period around a 
condition of "synchronism," for which the inclusion of higher order terms in Equation 50 would be 
necessary in reducing the drift data according to that theory. The data reduction in the section on 
reduction of orbits includes only those orbits falling within the minimum 66-day period around 
"synchronism." Further refinement of this reduction to include the 7 outside- of- synchronous 
orbits (according to the above criterion) will be made in the near future. This refinement is not 
expected to materially affect the results of this report. 

Summarizing, the results show that within reasonable excursions, the drift in longitude from 
a synchronous mean longitude is given approximately but explicitly in terms of time as 

/ \ A * 2 

K - k 0 - [A 22 sin 2 y 0 ) — * (54a) 

The corresponding mean daily radial drift is 


(A 22 ) 2 (At) 3 


1877 


(54) 


/ 1. 185 X 10~ 5 X 1877 \ 13 

V (23. 2) 2 x 10~ 12 / 


108 sid. days from 
synchronism. 


a 


l 



sin 2 y 0 
3m 


) 


At . 


(54b) 


Appendix E shows that Equations 54a and 54b are good approximations of the actual drift solutions 
if At is limited to the order of about 100 days for the probable order of magnitude of the earth’s 
dominant tesseral gravity harmonic. The study in Reference 4 comes to the same conclusion. It 
is clear that the two drift periods for Syncom Et, August to November 1963 and December to March 
1963-1964, in the close neighborhood of momentary synchronism at 54.8°W and 59.2°W fulfill in 
large measure the requirements for accuracy of the explicit drift equations 54a and 54b. The 
actual analysis of the data relied most heavily on 54a, since the mean daily change of the semimajor 
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axis was more difficult to extract from the limited data because of sun and moon perturbations 
which had large variable effects on orbits with epochs at different anomalies from the ascending 
node of the satellite. Note that Equation 54a, differentiated twice, yields 

K = y = - A 22 sin2y 0 , ( 54 c ) 

which is sufficiently accurate for the greater part of the two separate drift periods for Syncom n 
(see Appendix E). Thus, knowing k from the actual data for the two drift periods, Equation 54c can 
be solved uniquely for A 22 and y 0 , giving the magnitude and phase angle of earth equatorial el- 
lipticity assumed to have caused the accelerated drifts in the vicinity of the two separated "syn- 
chronous" longitudes. 


DETERMINATION OF EARTH-EQUATORIAL ELLIPTICITY FROM TWO OBSERVATIONS 
OF DRIFT ACCELERATION AT A GIVEN LONGITUDE SEPARATION 

While it is true that two evaluations of the near- synchronous drift accelerations at separated 
longitudes will determine the amplitude (a 22 , and thus J 22 ) of the perturbing sine function, the 
phase angle (y 0 , and thus k 22 j will have multiple solutions. However, since it is known from the 
section on solutions of drift equations that the drift acceleration is always in the direction of the 
nearest extension of the minor equatorial axis of the earth, the proper quadrant of the phase angle 
of equatorial ellipticity always can be resolved. 

Given two independent near -synchronous drifts (in the sense discussed previously), whose 
momentarily synchronous longitudes (y 0 ) l and (y 0 ) 2 are separated by V\. Let the two drift 
accelerations at these two "synchronous" configurations be (r 0 ), and ( > : „ ) 2 • The drift accelera- 
tions may be determined from drift-data reduction according to the theory of Equation 50. From 
Equation 33, 


(r'o)i ~ (^22)1 sin 2 {yo)i 


(55) 


and 


('>'0)2 (^22)2 s i n 2 ^("y 0 )i + > 


(56) 


since (r 0 ) 2 -(r 0 )i = 


(A-o) 2 ~ (x 0 )i • Expanding Equation 56 and dividing by 55 gives 


~(>o) 2" 

(^22 ) 1 

bo) J 

_(^22 ) 2 _ 


j^cos 2Vk + sin 2W_ cot 



(57) 
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Solving Equation 57 for ( y Q ) x , 


(*o)l 


1 

2 tan 


' fa ) 2 (^22)1 

(^o)l (^22)2 


1 2V\ 


; 2V\ 


( 58 ) 


The quadrant of [y 0 ) 1 is either the first or the fourth, because drift acceleration is always in 
the direction of the nearest longitude extension of the earth's minor equatorial axis.* Once the 
minor axis is located by Equation 58, the absolute value of J 22 in the earth's triaxial gravity field 
can be determined through Equations 55 and 30a, for example, as 


72t7 2 


(^22 )] 

[ R o/K) i ] 2 


(7o)i 


cos 2 i. + 1 


cos 2 i . + 1 


72tt 2 sin2(y 0 ), |^R 0/ /(a s ) jJ 2 3 


(59) 


Note that the units of (y 0 ) i in Equation 59 must be those of radians per sidereal day 2 , so that J 22 
will be dimensionless. Note also that in Equation 58, using the result of 30a, 


(^ 22)1 

(^22)2 


( a s ) 2 2 COs2 *1 + 1 
( a s) 1 cos 2 i 2 + 1 


(60) 


Using Equation 58, since (\ 0 )i is known from the data reduction (the geographic longitude of the 
" synchronous" configuration), the geographic longitude of the nearest minor axis location can be 
calculated as 


^22 (^ 0)1 (^ 0)1 ■ (60a) 

Similarly, the geographic longitude of the nearest major equatorial axis location can be calculated 
from 


^22 " (^0)1 " ( 7 0 )i + 90 degrees 


(60b) 


(see Figure 3). 

Following the theory of Wagner, f the difference in major and minor equatorial radii of the 
earth’s triaxial geoid a 0 -bQ is related to the gravity constant J 22 by 

a 0 ' b 0 - “ 6R 0 J 22 * (61) 

*Thus, the quadrant of the tan 1 in Eq» 58 can be most singly resolved by choosing the first or second quadrants if (^q) j < 0, and the third 
or fourth if (?o)l > 0* 

Wagner, op. cit. (See footnote, p. 7). 
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REDUCTION OF 27 SYNCOM II ORBITS TO DETERMINE THE EARTH’S 
EQUATORIAL ELLIPTICITY 


Appendix D tabulates the 27 Syncom II orbits from which the reduction below was made. 

Table 1 gives the estimated ascending equator crossings nearest to the epoch of these orbits. 

These were calculated by hand and therefore are listed only to 0.01 degree and 0.01 day. The 
technique used was to locate from the Nautical Almanac the geographic longitude of the ascending 
node at epoch through the reported right ascension of the ascending node for the orbit, and the hour- 
angle of the vernal equinox calculated at epoch. The geographic longitude of the ascending equator 
crossing was then estimated by turning the earth backward or forward through the orbit angle 
from the ascending node to the satellite at epoch. This latter quantity was estimated as + M, or 
360° - O + M), for the near-circular orbit of Syncom II. A correction factor to this orbit angle— 
the ratio of the satellite’s period to the earth’s sidereal period— was applied for orbits whose 
period was sufficiently different from the earth’s. The nodal longitude at epoch, plus or minus this 
reduced nodal excursion angle, is the estimated ’’ascending equator crossing nearest to epoch’’ re- 
ported in Table 1 (see Appendix D for an example of this calculation). 


Table 1 


Estimated Ascending Equator Crossings Nearest the Epoch of 27 Syncom II Orbits 


First Drift 
Orbits 

Time from 

Ascending Equator 

Second Drift 
Orbits 

Time from 

Ascending Equator 

20.0 Aug. 1963 

Crossing (degrees 

26.0 Nov. 1963 

Crossing (degrees 

(days) 

west of 50.0°W) 

(days) 

west of 50.0°W) 

1-1 

2.12 

4.89 

2-1 

1.86 

9.17 

1-2 

7.11 

4.83 

2-2 

7.84 

9.17 

1-3 

11.09 

4.78 

2-3 

13.83 

9.22 

1-4 

16.08 

4.74 

2-4 

20.81 

9.38 

1-5 

20.07 

4.77 

2-5 

41.75 

10.15 

1-6 

23.06 

4.78 

2-6 

44.74 

10.36 

1-7 

28.05 

4.85 

2-7 

55.71 

11.02 

1-8 

31.04 

4.90 

2-8 

64.69 

11.76 

1-9 

38.02 

5.06 

2-9 

71.67 

12.32 

1-10 

42.01 

5.09 

2-10 

76.66 

12.81 

1-11 

48.99 

5.45 

2-11 

83.64 

13.49 

1-12 

54.97 j 

5.68 




1-13 

62.95 

6.09 




1-14 

70.93 

6.60 




1-15 

77.91 

7.14 




1-16 

83.90 

7.61 





-100.0 

First free-drift 
period ends at an 
ascending equator 
crossing of =9.15° 
west of 50.0° West 
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Table 2 gives the Goddard- reported semimajor axes for these 27 orbits. Truncating Equa- 
tions 50 and 53a at their first right-hand terms gives 


AX. (longitude drift from “synchronism”) 


- (A 22 sin 2 y 0 ) 


(At) 2 
2 ’ 


(62) 


a 2 (semimajor axis change from “ synchronism”) 


sin2y 0 (At) 

A 22 3tt 


(63) 


Let the drift time be given from a certain arbitrary base time by T; let T 0 be the time of "syn- 
chronism" from the base time. Let the drift be given from a certain arbitrary geographic longi- 
tude by k and let k 0 be the geographic longitude from this base longitude of the "synchronous" 
configuration. Then 


At = T - T 0 , 


and 


A k = k - k Q . 


Table 2 


Goddard- Reported Semimajor Axes for 27 Syncom H Orbits 


First Drift 
Orbits 

Time from 
20.0 Aug. 1963 
(days) 

Semimajor Axis 
(42160.0 + data; km) 

Second Drift 
Orbits 

Time from 
26.0 Nov. 1963 
(days) 

Semimajor Axis 
(42160.0 + data; km) 

1-1 

2.27 

4.58 

2-1 

2.04 

5.89 

1-2 

6.71 

4.52 

2-2 

8.00 

7.20 

1-3 

11.00 

6.02 

2-3 

14.00 

7.18 

1-4 

16.00 

6.39 

2-4 

20.71 

8.17 

1-5 

20.00 

6.35 

2-5 

41.71 

8.01 

1-6 

23.08 

6.55 

2-6 

44.25 

9.90 

1-7 

28.08 

6.70 

2-7 

55.88 

11.43 

1-8 

31.08 

7.42 

2-8 

64.83 

11.91 

1-9 

38.08 

7.51 

2-9 

71.67 

12.89 

1-10 

42.08 

8.88 

2-10 

76.79 

13.31 

1-11 

49.08 

9.14 

2-11 

83.71 

14.89 

1-12 

55.08 

9.78 




1-13 

63.08 

11.51 




1-14 | 

71.00 

11.09 




1-15 

78.00 

12.15 




1-16 ! 

84.21 

12.51 





= 100.0 

First free drift 
period ends with a 
semimajor axis of 
=42174.5 km 
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With these changes, Equations 62 and 63 become [noting that a, = (a - a s ) / /a s 


A 22 sin 2y 0 

h = h 0 2 


T 2 - 2TT 0 + T 0 2 


or 


T o ( A 22 sin 


+ T (T 0 A 22 sin 2 y 0 ) + T 2 (- 


A 22 sin 2y 0 \ 


; . 


(64) 


(Note that, from Equation 64, \ = - A 22 sin 2y 0 = y 0 from Equation 33. This result is valid only for 
orbits sufficiently close to "synchronous," as discussed previously.) 


a - + a 


(i^) (T . TJ 


s s 


or 


/ T 0 a 22 sin 2y 0 \ /a s A 22 sin 2y, 

Y 35? / + T \ 


0 


(65) 


Equations 64 and 65 may be written with determinable coefficients as 


K = d 0 + d, T + d 2 T 2 


e 0 + e t T 


( 66 ) 

(67) 


where 


dp ho 


T 0 2 A 22 sin 2 r 0 


dj + A 22 Tq sin 2y 0 


A 22 sin 2y 0 


(>- 


T q A 22 sin 2y 0 ^ 
3?7 / 


e l a s A 22 


sin 2 y 0 

3?7 




V (68) 
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From Equation 68, 


T o - - dj/iWj , (69) 

( d ,) 2 

“0 = d 0 - Tap ’ (70) 

r'o r - A 22 sin 2y 0 = 2d 2 • (71) 

Alternately, and as an internal check on the theory of the coupling of the drift and orbit expansion, 

- 3me l 

r’o = ” A 22 sin 2 ?0 " i ’ 


implying 


d 2 - fce^a s . 


(72) 


In Equation 72, the units of d 2 must be radians per sidereal day 2 , and the units of e t must be 
length per sidereal day so that the equation will be dimensionally correct. The semimajor axis at 
the "synchronous" configuration is calculated from Equation 67 for T = T 0 : 


a s - e o + e i T o • 


(73) 


For the first drift period (orbits 1-1 through 1-16), the best estimates (in the "least squares" 
sense) of the coefficients (d)j and (e), , obtained by fitting Equations 66 and 67 to the data in Tables 1 


and 2, have been found to be 


( d o)i 

= 4.941 ± 0.018 degrees. 

( d l)l 

= - 0.0216 ± 0.0010 degree/solar day, 

( d 2 ). 

= (6.37 ±0.11) x 10~ 4 degree/solar day 2 
= (6.33 + 0.11) x 10‘ 4 degree/sid. day 2 , 

( e o)l 

= 4.35 ± 0. 19 km, 

( e l)l 

= 0.0993 + 0.0042 km/solar day 

= 0.0990 ± 0.0042 km/sid. day. 
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The mean value of the inclination during this period was 


ij = 33.018 ± 0.005 degrees. 


From Equation 69, 


( T o)l 


16.95 _ 


1.09\ 
1. 05y 


days from 20.0 August 1963. 


From Equation 70, 


( \ 0 ) i = 4.76 ± 0.03 degrees west of 50.0°W long. 


From Equation 71, 

( 7 * 0)1 " " (1.27 ±0.02) X 10 " 3 degree/solar day 2 = - (2.20 ±0.04) x 10 “ 5 rad/sid. day 2 . 

From Equation 73, and the above value of (T 0 )! , 


(a g ) j = 42166.0 ±0.2 km. 

For the second drift period (orbits 2-1 through 2-11), the best estimates (in the ’’least squares’ 1 
sense) of the coefficients d 2 and e 2 , obtained by fitting Equations 66 and 67 to the data in Tables 1 
and 2, have been found to be 

(d 0 ) 2 “ 9.156 ± 0.017 degrees, 

(di ) 2 = - (0.0030 ± 0.0010) degree/solar day, 

(d 2 ) 2 " ( 6 . 59 ±0.11) x 10 " 4 degree/solar day 2 , 

= ( 6.55 ± 0 . 11 ) x 10~ 4 degree/sid. day 2 , 

(e 0 ) 2 = 5.70 ± 0.42 km, 

(e 1 ) 2 - 0.0994 ± 0.0080 km/solar day 

~ 0.0990 ± 0.0080 km./sid. day. 

The mean value of the inclination during this period was 

i 2 = 32.851 ± 0.010 degrees. 
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From Equation 69, 


(To) 2 - 2.3 ± 0.8 days from 26.0 November 1963. 


From Equation 70, 


(^ 0)2 ~ 9-15 ± 0.02 degrees west of 50.0°W long. 

From Equation 71, 

(> 0)2 ~ ” (1.32 ±0.02) x 10 -3 degree/solar day 2 = - (2. 29 ±0.04) x 10~ 5 rad/sid. day 2 . 

From Equation 73, and the above value of (t 0 ) 2 , 


( a *± 


42165.9 + 0.4 km. 


(See Figure 5 for a graph of this orbit data and reduction for the two drift periods.) Combining 
the above results of the two free-drift periods, from Equation 60, 


( A 2 2 ) 

(M 


= (42165.9 + 0.4/42166.0 ±0.2) 2 


cos 2 (33.018± .005) + 1 
cos 2 (32.851 ± .010) + 1 


= 0.99845 ± .00014. 


The longitude separation between the two drift periods is given by 


V\ = (\ 0 ) 2 - (\ 0 ), = [-(59. 15 ±0.02)] - [-(54.76 + 0.03)] degrees 


= - (4.39 ±0.05) degrees geographic longitude. 


Thus, 


2V\ = - (8.78 ±0.10) degrees geographic longitude. 


Therefore, from Equation 58, the location of the minor equatorial axis with respect to the "syn- 
chronous" longitude during the first free-drift period (54.76 ± 0.03 degrees west of Greenwich) is 


( y °)i 


2 tan“ 


sin [- (8.78+0. 10)1 


= (»-J) 


1.27 + a0 2 (0.99845 +0.00014) - cos [- (8.78 ± 0. 10)] 
degrees east of the minor equatorial axis. 
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Figure 5— Drift of the ascending equator crossing and growth of the semimajor axis of Syncom il 



I 


From Equation 60b, the best estimate of the location of the major equatorial axis is 


"22 


55 


/ + 5\ / +5\ 

154 _ g I + 90 - - I degrees geographic longitude. 


( 74 ) 


From Equation 59, the best estimate of the triaxial gravity coefficient J 22 is 


» 22 


72rr 2 sin 2 


(«•!?) 


( 2. 20 ± 0. 04) x 10 


rs 


(6378.2/42166.0 ±0.3) 


2 cos 2 (33.018 ±0.005) ± 1 


/ + 0.07\ 

= " \± 67 -O.03J * 10 • (75) 

The mean equatorial radius, taken as R 0 = 6378.2 km, is a compromise for a number of cur- 
rently used values. It is stated above without error. The likely error in (a s ) 2 has been increased 
arbitrarily by 0. 1 km to account for the likely uncertainty in R 0 . 

The results of the simulated Syncom n trajectories for these drift periods (Appendix C)— with 
equatorial ellipticity parameters close to Equations 74 and 75 included in the particle program of 
the simulation— shows that sun, moon, and earth zonal gravity over the 3 -month drift periods 
studied does not substantially affect the simple theory of this reduction. 

However, the simulations show small biases in the simple reduction for J 22 and \ 22 from 
Equations 62 and 63 on the basis of an elliptical earth equator only. 

The final reduced geodetic parameters adjusted for cumulative sun, moon, and earth zonal 
gravity effects are: 

J 22 = - (1.70 ±0.05) x 10' 6 ( 76a ) 

(where J 22 is adjusted for all but higher order longitude earth-gravity effects), and 


x 22 = - (19±6)° (7 6b ) 

(where \ 22 is similarly adjusted; see Appendix C). As noted, higher order longitude earth-gravity 
effects have not been taken account of in these final results. 

As can be seen from Figure Al, (the earth's longitude gravity field at the equator according to 
a recent estimate) the effect of higher order gravity on the results of Equation 76 is not expected 
to exceed 15 percent on J 22 and a few degrees on ± 22 . Figure 6 shows a comparison of the 24- 
hour altitude longitude gravity field along the equator as derived herein from Syncom n drift over 
Brazil and three other recent estimates from lower altitude satellite and surface gravity data. 
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0 90 180/- 180 -90 0 


GEOGRAPHIC LONGITUDE (degrees) 


Figure 6— The earth's longitude gravity field at "synchronous altitudes" on the 
equator, according to recent estimates on different bases (see Table Al). 


Using the above estimate of J 22 from observations on Syncom II drift, the difference between 
the major and minor equatorial radii of the triaxial geoid is, by Equation 61, 


a o - b 0 - 65 ± 2 meters - 213 + 6 ft. 

Comparing the deviation due to earth ellipticity with other higher order earth- gravity devia- 
tions (Appendix A*), we note that the above figure implies a maximum deviation from the mean 
earth sphere, due to the ellipticity of the equator, of 

AR 0 = 107 + 3 ft. 


STATION-KEEPING REQUIREMENTS FOR EQUATORIAL 24-HOUR SATELLITES 

Following the calculation first made by Fleig, 1 maximum station- keeping requirements for a 
synchronous equatorial satellite in a purely triaxial earth field occur at a longitude 45 degrees 


♦Also, Wagner, op. c it. (See footnote, p. 7). 

h*'leig, A., "Effect of Nonsphericity of Earth’s Gravitational Potential on an Equatorial Synchronous Satellite,” NASA Goddard Space 
Flight Center, Guidance and Control Section (Code 622) Report 43, June 15, 1962. 
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from the major axis of the elliptical equator. The triaxial longitude perturbation force which 
must be continually counteracted by on-board propulsion devices to achieve perfect station keep- 
ing is given in Appendix A (Equation A4) as 


( ) 2 2 



sin 


2 



(77) 


The maximum perturbation force occurs at longitudes where \-\ 22 = ±45, ±135 degrees; or, for 
the reported adjusted value, ^ 22 = -19 degrees; maximum perturbation forces occur at 


| K - 64, - 154, 116, 26 degrees . 

| 

| For the reported adjusted value; J 22 - -(1.7) x 10“ 6 , with 

Me/ ( a s) 2 = 0.735 ft/sec 2 . (R 0/ a s ) 2 - 0.0229 , 

the maximum yearly impulse per unit mass, or velocity requirement for perfect station keeping 
(at these longitudes), is: 

(F^) 22 m3 x x T ( 1 y ear ) = , max = 0 . 735 x 0 . 0229 x 86 , 400 X 365 X 6 x 1 . 70 x 10" 6 

“5.36 ft/sec/yr . (79) 


The result in Equation 79 is strictly true only for a triaxial earth field. As Reference 4 shows, 
and the combined geoid of Kaula (1964) confirms, the absolute maximum longitude perturbation 
(at 24-hour altitudes) in a full earth field may increase by as much as 15 percent over that in the 
simple triaxial field. According to the combined geoid of Kaula (1964), this absolute maximum 
longitude perturbation occurs over Indonesia (see Figure Al, Appendix A). 

For synchronous satellites to be kept permanently in this quadrant may require, conservatively, 


£V t . - 6.17 ft/sec/yr . (80) 

However, this latter result may not be conservative enough if there is an additional bias in the 
reported value of J 22 (in Equation 76) due to the neglect, in that reduction, of these same higher 
order effects. A study of Kaula's geoid of Figure Al, Appendix A, certainly indicates that such a 
bias would be present for a J 22 -only reduction from the full field effects at the two longitudes 
-54.8 and -59.2 degrees. An investigation is underway by the author of these probable higher order 
effects on the reported Syncom II geodetic reductions from drift over Brazil (Equation 76). 


DISCUSSION OF RESULTS 


Some discussion is in order on the validity of the reported geodetic parameters (Equa- 
tion 76). These adjusted parameters, without accounting for higher order earth-gravity 
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effects, show a somewhat stronger and west- shifted equatorial ellipticity than the most 
recent satellite geoids (see Table Al, Appendix A). Around the equator, the reported Syn- 
com II determination yields a longitude gravity field at 24 -hour altitudes in reasonable agreement 
with recent satellite and surface gravimeter determinations (see Figure 6). The basic assumption 
behind the reported Syncom n determination is that during the analyzed drift periods the satellite 
moved free of all external influences except sun, moon, and earth gravity. 

It seems very unlikely that other small influences on Syncom II motion, such as (1) micro- 
meteorite collisions, (2) possible leaking gas from the satellite, (3) solar radiation and solar wind 
pressure, (4) magnetic field interactions, and (5) planetary gravity, could have had appreciable 
selective influence on the observed accelerated drift of the satellite. The fact is that the simulated 
Syncom n drifts (Appendix C) in a sun, moon, and earth-gravity field only reproduce the observed 
drifts in all particulars to a degree which seems to preclude other influences except small random 
ones. 

From injection to mid- August 1963, there appeared to be some question as to possible leaking 
of gas on board Syncom II. Monitored gas pressure levels in some of the propellant tanks de- 
clined erratically between correction maneuvers in this period. However, levels were reasonably 
constant in the analyzed long free-drift periods. The evidence of erratic and declining gas pres- 
sure levels on board Syncom n is thought to be due, at least in part, to pressure transducer 
calibration errors. 

Sun, moon, and zonal earth-gravity influences on Syncom H T s long-term drift acceleration, as 
well as model errors made in the assumptions of orbit circularity and unchanging inclination, have 
been accounted for in the simulated drift reductions. As yet, only two sustained longitude samples 
of the earth ! s gravity field have been taken by a 24-hour satellite. It is noted from the evidence 
of Figure Al (Appendix A) and Reference 4* that the ”true value” of the magnitude of equatorial 
ellipticity, when all higher order earth-gravity effects are accounted for, may be as much as 15 
percent higher than that reported herein. With the next few samplings at reasonably separated 
longitudes, a considerable refinement in the estimate of equatorial ellipticity (as affected by higher 
order gravity) should be possible. Along with this refinement should come an order-of- magnitude 
estimate of third- (and higher) order longitude earth gravity. Much geodetic information remains 
to be revealed by closely following the drift of Syncom II and other 24-hour satellites soon to be 
launched by NASA. 

CONCLUSIONS 

From this investigation into longitude-gravity-caused Syncom II drift, the following conclu- 
sions are made: 

1. Longitude- dependent earth gravity exists and must be reckoned with in the long-term opera- 
tion of 24 -hour satellites. 


•Also: Private communication from W. M. Kaula, 1964, ^Theory of Satellite Geodesy,” (unpublished text for a course at the Institute of 
Geophysics and Planetary Physics, UCLA, Los Angeles, Calif.) 
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2. As far as 24-hour satellites are concerned, the earth can be assumed to be only a triaxial 
ellipsoid with the following measures of equatorial ellipticity: 

(a) J 2 j = -(1.70 +0.05)* 10" 5 6 , 

corresponding to a 65 ±2 meter difference between major and minor equatorial radii; 

(b) = - (19 ±6)° , 

locating the geographic longitude of the major equatorial axis. 

3. Higher order earth- gravity effects appear to distort the simple triaxial field at 24- hour 
altitudes, by about 10 to 15 percent, at most, around the equator. A few of these effects of the 
third- and fourth-order should be well discriminated by long-term observations of future 24-hour 
satellites. 

4. With negligible higher order earth effects assumed, the maximum longitude gravity station- 
keeping requirements for the 24-hour satellite are 5.36 ft/sec/yr. 
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Appendix A 

Earth-Gravity Potential and Force Field Used in this Report: 

Comparison with Previous Investigations 

The gravity potential used as the basis for the data reduction in this study is the exterior po- 
tential of the earth derived* for geocentric spherical coordinates referenced to the earth's spin 
axis and its center of mass. The infinite series of spherical harmonics is truncated after J 44 . 

The nontesseral harmonic constants J 20 , j 30 , and j 40 are derived from Reference 8. 

The earth radius R 0 used in this study is 

R 0 = 6378.388 km. 

The earth's gaussian- gravity constant used is 

M E = 3. 9862677 x 10 s km 3 /sec 2 . 

Neither of these values, taken from Reference 9, nor the "zonal geoid" of Reference 8 is felt to be 
the most accurate known to date. They are the values used by the GSFC Tracking and Data Sys- 
tems Directorate to calculate the orbit elements of Syncom n from radar and Minitrack observa- 
tions. They were chosen to insure consistency between the data of this study and these published 
orbits, inasmuch as the "triaxial" reduction for which this study has been undertaken is not sig- 
nificantly affected by the probable errors in these values. The second-order tesseral harmonic 
constants used in the simulation studies were 

J 22 = - 1.68xl0' 6 , 

/V 22 = - 18 degrees . 

These are the values shown on the "tesseral geoid" below (for the j 22 harmonic). At a later point 
in the analysis, the slightly different values reported in the abstract were estimated. The most 
accurate "zonal geoid" is probably that of Kozai(see Referenced), with the following earth constants: 
R 0 = 6378.2 km, 

J 20 = 1082.48 x 10' 6 , J 3() = - 2.56 x 10~ 6 , J 40 = - 1.84 x 10~ 6 , 

M e = 3.98603 x 10 s km 3 /sec 2 . 

The earth's gravity potential (to fourth order, probably sufficient to account for all significant 
longitude perturbations on a 24-hour satellite) is given by Equation Al, which may be illustrated 
as follows (following Reference 4, Appendix B): 


♦Wagner, op, cit. (See footnote, p. 7.) 


35 






The earth-gravity field (per unit test mass), whose potential is Equation Al, is given as the gra- 
dient of Al, or 


F = r F r + \F X + = VV £ 


+ X 

r <9 r + r co s (p r9 X 


0^E 

r d(p 


(A2) 


or 

F r = {" 1 + ( R o/ r ) 2 [3/2J 20 (3 sin 2 <t>- l) + 9J 22 cos 2 cpcos 2 (\ - A. 22 ) 

+ 2(R 0/ /r) J 30 ($ sin 2 <£ - 3) (sin <£) +6(R 0 /r) J 31 '5 sin 2 ^ “ l) cos^cos (\- k 31 ) 

+ 60(R 0 /r) J 32 cos 2 <p sill <£ cos ~ ^ 32 ) + 60 (R 0 /r) J 33 cos 3 <p cos 3 - \ 33 ) 

+ 5/8(R 0 /r) 2 J 40 (35 sin 4 - 30 sin 2 0 + 3 ) 

T 25/2 (R 0 / r ) 2 J 41 (7 sin 2 0 - 3 ) cos 0 sin 0 cos (x ~ X 41 ) 

+ 75/2(R 0 /r) 2 J 42 (7 sin 2 0- l) cos 2 0 cos 2 (x - X 42 ) 

+ 525 (R 0 /r) 2 J 4 3 cos 3 <p sin 0 cos 3 (X - X 43 ) + 525(R Q /r) 2 J 44 cos 4 <0 cos 4 (^M}’ (A3) 


0 E . 

F A = — (R 0 /r) 2 |6J 22 cos<£ sin 2 (\ ~\ 22 ) + 3/2 (R 0 /r) J 3l (5 sin 2 - l) sin(\-\ 31 ) 

+ 30(R o/ r) J 32 cos</; s i n 0 sin 2 ( X “ ^ 32 ) +45 ( R o/ r ) J 33 cos 2 0sin3(A.-\ 33 ) 

+ 5/2 (R 0 /r) 2 J 41 (7 sin 2 4> - 3) sin sin(\ - \ 41 ) + 15 (R 0 /r) 2 J 42 (7 sin 2 0 - l) cos <p sin 2 (\ - \ 42 ) 
+ 315(R 0 /r) 2 J 43 cos 2 0 s i n 0 s i n 3 (a. - X 4 

+ 420 (R 0 /r) 2 J 44 cos 3 0 sin 4 (\ - \ 44 )j , ( A4 ) 


Me r 

= — (F 0 /r) 2 |“ 3J 20 sin i cos d-r 6J 22 cos 0sin <?cos 2(>v- X 22 ) 

- 3/2(R 0 /r) J 3Q (5 sin 2 <p- l) cos P + 3/'2(R 0 /r) J 3l (15 sin 2 - ll) sin <£cos('. - \ 31 ) 

+ 15 (R 0 /r] J 32 (3 sin 2 <£- l) cos (pc os 2 (k~ \ 32 ) 

+ 45(R 0 /r) J 33 cos 2 <p sin <p cos 3 (\ - ^ 33 ) - 5/2(R 0 /r) 2 J 40 (7 sin 2 0 - 3 ) sin<± cos p 

+ 5/2(R 0 /r) 2 J 41 (28 sin 4 (+-27 sin 2 p + 3 ) cos(^ ~ \ 41 ) 

+ 30(R 0 /rJ 2 J 42 (7 sin 2 <p~ 4) cos <£sin pcos 2{\ - 

+ 105 (R 0/ r ) 2 J 43 (4 sin 2 l) cos 2 pc os 3(/v~ \ 43 ) 

+ 420(R 0 /r) 2 J 44 cos 3 <£sinc cos4(\“\ 44 )j (A5) 

The actual sea-level surface of the earth is to be conceptualized through Equation A1 as a sphere 
of radius 6378 km, around which are superimposed the sum of the separate spherical harmonic 
deviations illustrated. To these static gravity deviations, of course, must be added a centrifugal 
earth- rotation potential at the earth's surface, to get the true sea-level surface.* 

Table A1 gives tesseral coefficients for this earth- gravity field form as reported by geodesists 
from 1859 to 1964. Figure A1 presents the longitude gravity field at the altitude of the 24-hour 
satellite, around the equator, according to the geoid of Kaula- Combined (1964). Figures A2 through 
A8 give the earth surface contours of seven recent geoids whose individual tesseral coefficients 
are listed in Table Al. 


* Wagner, op, cit, (See footnote, p. 7.) 
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Table A1 (Continued) 
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LONGITUDE PERTURBATION 
IN FULL EARTH FIELD TO 4TH ORDER 
OF KAULA- COMBINED (1964) 



GEOGRAPHIC LONGITUDE (degrees) 













Appendix B 


Expressions for the inclination Factor 


Equation 25 in the text gives the inclination factor in the drift causing tangential perturbation 
(due to equatorial ellipticity) on a 24- hour satellite with a near- circular orbit as 

sin2 1 

F(i) = cos i + 2 ; (Bl) 

is the absolute value of the maximum longitude excursion of the figure- 8 ground track of the 
24-hour satellite (with a near-circular orbit) from the geographic longitude of the nodes. 

From Equation 18, this longitude excursion function is 

A\. = K - K g = tan " 1 (cos i tan#) - 6 . (B2) 

Differentiating Equation B2 with respect to the argument angle 8 , the minimax excursion 
arguments are found from 


d( A\) 
d 6 


0 


cos i sec 2 6 
1 + cos 2 i tan 2 0 


- X . 


Solving Equation B3 for Sin 6 at 


sin %A. 


min i max 


( cos i + 1) 1/2 


(B3) 


from which 


tan 6 


AX ,m in irnax 


- (sec i) 


1/2 


(B4) 


Equation B4 in B2 gives 


A\. = tan 1 (cos i) I/2 - tan" 1 (sec i) 1/2 . 

mini max v ' v ' 

Thus, since only the absolute value of AX roinimax is required, 

A\ max = tan' 1 (sec i) 1/2 - tan -1 (cos i) 1/2 , ^ (B5) 
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where the tan 1 is to be taken in the first quadrant. For example, for i = 30 degrees, Equation 
B5 evaluates the maximum excursion as 

A\ - 47.0-42.9 - 4.1 degrees . 

max 13 

The nodal argument angle at this maximum longitude excursion is 


#(at AX. max ) - ±42.9 degrees from the nodes . 


The assumption in Equation 20 (in text) that the excursion in longitude from the ascending node 
could be approximated by 


A\ = -A\ sin 26 

max 


predicts the maximum excursion argument as 


6 (at A\ max ) - ±45 degrees from the nodes . 

This discrepancy in the assumed longitude excursion function is not serious until i > 45 degrees, 
as simulated trajectories with variable inclination have borne out. 

Equation B5 can be written as 

A\ max + tan' 1 (cos i) 1/2 r - tan -1 (sec i) 1/2 , 


from which 


. x A\ + (cos i) 1/2 

tan[A\ + tan -1 (cos i) 1/2 J = (sec i) 1/2 - * ; r~rr- , 

\ max V J 1 K } 1 - AX (cos i) 1/2 

max' ' 


for i < 45 degrees. Solving Equation B6 for AX max , 

^ 1 - cos i A 1 - cos i 

max 2(cos i ) 1 7 2 1 + cos 1 


for i small. 

Thus, the inclination factor becomes approximately 


F( i ) - cos i + 


sin 2 i ( 1 - cos i ) 
2( 1 + cos i ) 


cos 1 + 


( 1 ~ cos i) 2 


2 cos i + 1 - 2 cos i + cos 2 i 


cos 2 i + 1 
2 


(B6) 


(B7) 


(B8) 
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For example, for i = 30 degrees. 


F(i) from Eq. B1 " 0.86603 + 4.1/8x 57.3 = 0.8750, 

F(i) from Eq . Be = 0.8750. 

The agreement of F( i ) from forms B1 or B8 is good to the third decimal place as long as the in- 
clination is less than 45 degrees. At inclinations higher than 45 degrees, however, the drift theory 
following Equation 20 (in text) begins to break down because ^ max is no longer a small angle. 
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Appendix C 


Reduction of Simulated Particle Trajectories for 
Earth-Equatorial Ellipticity 

Tables Cl and C2 present data taken from two numerically integrated particle trajectories of a 
triaxial earth in the presence of the sun and moon's gravity field. Only perturbed equations of 
motion from a periodically rectified Keplerian reference orbit are actually integrated by the 
digital computer program (called ITEM at Goddard Space Flight Center). For the three months' 
real orbit time of these trajectories, the accumulated truncation and roundoff error in the nu- 
merical integration is believed to be negligible for the purposes of this reduction. The initial posi- 
tion and velocity conditions for these simulated trajectories were the same as those reported for 
the "actual" Syncom n orbits 1-2 (for the trajectory of Table Cl) and 2-3 (for the trajectory of 
Table C2). The program used the earth gaussian- gravity constant 

f ij, = 3.9862677 x 10 s km 3 /sec 2 , 

which is the gravity constant used by the GSFC Data and Tracking Systems Directorate in comput- 
ing the elements of satellite orbits from radar and Minitrack observations. The best estimates 

Table Cl 

Data from Simulated Trajectory Beginning with the Elements of Syncom n Orbit 1-2 


(J 22 = -1-68 x 10~ 6 , R 0 = 6378.388 km,y 22 = -108.0: input into trajectory program) 


Time from 
26.709 Aug. 1963 
(solar days) 

Ascending Equator Crossing 
(degrees west of 50.0°W 
geog. long.) 

Semimajor Axis 
(42160.0 + data; km) 

Inclination 

(32.0 + data; degrees) 

2.390 

4.816 

5.27 

1.089 

8.374 

4.783 

7.09 

1.072 

14.358 

4.792 

6.01 

1.056 

20.341 

4.861 

8.12 

1.043 

26.324 

4.954 

7.13 

1.025 

32.308 

5.101 

8.98 

1.019 

38.292 

5.291 

8.31 

0.997 

44.276 

5.537 

9.67 

0.991 

47.268 

5.678 

10.38 

0.983 

50.260 

5.821 

10.03 

0.972 

53.253 

5.975 

9.42 

0.967 

56.245 

6.144 

9.74 

0.966 

59.237 

6.326 

11.09 

0.960 

62.229 

6.522 

11.94 

0.957 
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Table C2 

Data from Simulated Trajectory Beginning with the Elements of Syncom II Orbit 2-3 


( J 22 = -1.68 x 10 6 , R 0 = 6378.388 km,y 22 = -108.0: input into trajectory program) 


Time from 
10.0 Dec. 1963 
(days) 

Ascending Equator Crossing 
(degrees west of 50.0°W 
geog. long.) 

Semimajor Axis ' 

(42160.0 + data; km) j 

Inclination 

(32.0 + data; degrees) 

0.823 

9.243 

6.88 

0.881 

5.809 

9.351 

7.30 

0.881 

10.796 

9.495 

9.41 

0.877 

15.783 

9.666 

8.15 

0.864 

20.769 

9.885 

9.85 

0.864 

25.756 

10.134 

10.60 

0.850 

30.743 

10.401 

9.95 

0.842 

35.730 

10.708 

11.81 

0.841 

40.717 

11.044 

12.18 

0.825 

45.704 

11.412 

11.58 

0.816 

50.692 

11.830 

13.93 

0.808 

55.679 

12.259 

13.11 

0.790 

58.672 

12.534 

13.07 

0.785 

60.667 

12.724 

13.69 

0.784 


(in the n least squares" sense) of the coefficients (d) sl and (e) sI , obtained by fitting the drift and 
orbit expansion text equations 66 and 67 to the data in Table Cl, have been found to be 


K) ,i 

( d i).i 

( d 2 )s> 

( e °) sl 

C e 0sl 


a 4.84110.004 degrees , 
n “ ( 1 . 22 ± 0. 03) x 10 -2 degree/solar day , 

n (6. 303 ± 0. 038) x 10 -4 degree/ solar day 2 , 

- ( 6. 268 ± 0. 038) x 10" 4 degree/sid. day 2 , 
n 5.45 ±41 km , 

- 0.091 ±0.010 km/ solar day , 
n 0.091 ±0.010 km/sid. day . 


The mean value of the inclination during this first simulated trajectory period was 


i sl - 33.005 ±0.003 degrees . 
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From Equation 69, 

( T o)sl = 9-68 + 0.30 days from 26.709 August 1963 . 

From Equation 70, in the text, 

, , / + ,008\ 

(— o)st ~ (4.782 - . 007 J degrees west of 50°W long . 

From Equation 71, 

( y q ) s j = ~ ( 1 . 261 ± 0 . 008) x 10' 3 degree/sol ar day 2 
= - (2. 188 ± 0.013) x 10 ' 5 rad/sid. day 2 . 

From Equation 73, and the above value of (T 0 ) sl , 

( a s)si = 42166.3 + 0.4 km . 

The best estimates (in the "least squares" sense) of the coefficients (d ) s2 and (e) j2 , obtained by fit- 
ting the drift and orbit expansion Equations 66 and 67 to the data in Table A 2, have been found to 

be 

{d 0 ) S 2 = 9. 224 ± 0.004 degrees 

(d j ) s 2 = (1.830 ± 0.028) x 10' 2 degree/solar day 

(d 2 ) s2 = ( 6 . 501 ±0.042) x 10 ' 4 degree/ solar day 2 

= (6.465 ±0.042) x 10 ' 4 degree/sid. day 2 

( e o )*2 = 7. 19 ±0.37 km 

( e l) S 2 = o.lll ±0.010 knv/ solar day 
= 0.111 ± 0.010 km/sid. day . 

The mean value of the inclination during the second simulated trajectory period is 

i s2 = 32.836 ±0.003 degrees . 
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From Equation 69, 


(T 0 ) s2 - - ( 14. 07 ± 0. 30) days from 10.0 December 1963 . 

From Equation 70, 

. . / + . 009 \ 

(-o) s2 “ (9 095 - 009 / degrees west of 50.0 W long . 

From Equation 71, 

(■y'o) s2 = -( 1. 300 ± 0. 008) x 10~ 3 degrees/ solar day 2 

= - (2. 257 ±0.015) x 10' 5 rad/sid. day 2 . 

From Equation 73, and the above value of (t o ) s2 , 

( a s) s2 = 42165.6 ±0.5 km . 


A graph of these trajectory simulations is seen in Figure Cl. 

Combining the above results for the two simulated trajectories, from Equation 60, 


22) si 
(^ 22 ) s2 


(42165.6 ±0.5/42166.3 ±0.4) 2 


cos 2 (33.005 ±0.003) + 1 
cos 2 (32.836 ±0.003) + 1 


0.99840 ±0.00007 


V\ = (\)s2 ~ ( K o) si = - (59. 095 ± 0.009) - [- (54. 782 ±0.008)] , 

•••2V\ = - (8. 626 ± 0. 034) degrees geographic long . 


Therefore, from Equation 58, the location of the minor equatorial axis with respect to the "syn- 
chronous longitude" during the first simulated trajectory (54.782 ± 0.008 degrees west of Green- 
wich) is 

sin [- (8.626 ±0.034)] 

1.300 ± 0.008 r i 

r 261 ±0.008 (°- 99840 - 0. 00007) - cos [- (8. 626 ± 0. 034)] 




52.5 4 2.5 degrees east of the minor equatorial axis . 
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Figure Cl— Drift of the ascending equator crossing and growth of the semimajor axis of simulated Syncom li trajectories. 




From Equation 60a, the best estimate of the geographic location of the nearest extension of the 
equatorial minor axis from the simulated trajectory data is 

(7 22 ) s ~ ~ 54.8 - (52. 5 ±2. 5) = -(107.3 ±2.5) degrees geographic long . 

This compares well with the input value of (y 22 ) s = -108.0 degrees used to compute the simulated 
trajectories. From Equation 59, the best estimate of the triaxial gravity coefficient j 22 from the 
simulated data (according to the theory of this report) is 

, . _ ~(2. 188 ±0.013) x 1(T 5 

72 ir 2 [sin 2(52. 5 ±2.5)] (6378. 388/42166. 3 + 0. 4) 2 |^ P S . 2 .. ( 33 - 005 j °~ 003) + 

= - (1.64 ±0.03)x 10~ 6 . 

The mean equatorial radius used in the simulation is R 0 = 6378.388 km, the same used to com- 
pute the ’’actual" Syncom II orbits from the radar and Minitrack observations. 

The above value of(j 22 ) s compares reasonably well with the input value of (j 22 ) s = -1.68x 10~ 6 
used to compute the simulated trajectories. 

The model error implicit in the difference between the reduced and inputed geodetic coeffi- 
cients for the simulated trajectories warrants an adjustment of the J 22 and k 22 reported in the 
section on the reduction of the "actual" Syncom II orbits. The values below appear sufficient to 
cover all the uncertainties of this reduction for a triaxial earth: 

J 22 (reduced for a triaxial earth) - - ( 1. 70 ± 0. 05) x 10~ 6 , 

X- 22 (reduced for a triaxial earth) = -( 19 ± 6) degrees geographic long . 

As the introduction indicates, the principal effect unaccounted for in the reduction is the possible 
influence of higher order earth tesseral anomalies on the drift of Syncom n. Figure A1 indicates 
that the accumulated influence on synchronous satellites of all higher order anomalies is small 
compared with the second-order anomaly according to recent satellite geoids. Close and contin- 
uing observations on the drift of 24- hour satellites will clarify this assumption, and should be 
rewarded in time by revelation of the "tesseral" anomalies through third-order with an absolute 
precision almost as good as that reported here for the second-order effect. 
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Appendix D 


Basic Orbit Data Used in this Report 


The Orbit elements for Syncom n in Table D1 were calculated at the GSFC Tracking and Data 
Systems Directorate from radar and Minitrack observations made during the slew-drift periods 
from mid- August 1963 to February 1964. 

As an example of the estimation of the ascending equator crossing nearest to epoch, consider 
the orbit geometry at epoch (Figure Dl): 6 January 1964 at 17.0 hours Universal Time (orbit 2-5). 


Table Dl 


Syncom II Orbital Elements, August 1963 to February 1964 


Orbit 

Epoch (Universal Time) 
(year-month-day- 
hour-min-sec) 

Semimajor 

Axis 

(km) 

Eccentricity 

Inclination 

(degrees) 

Mean 

Anomaly 

(degrees) 

Argument 

of 

Perigee 

(degrees) 

Right Ascen- 
sion of the 
Ascending 
Node 
(degrees) 

1-1 

63-8-22-6-1 2-8 

42164.58 

0.00018 

33.083 

24.126 

26.285 

317.569 

1-2 

63-8-26-17-0 

42164.52 

0.00016 

33.090 

190.841 

26.099 

317.454 

1-3 

63-8-31-0-0 

42166.02 

0.00018 

33.062 

296.125 

30.073 

317.475 

1-4 

63-9-5-0-0 

42166.39 

0.00012 

33.064 

333.521 

357.756 

317.362 

1-5 

63-9-9-0-0 

42166.35 

0.00015 

33.048 

326.207 

9.077 

317.272 

1-6 

63-9-12-2-0 

42166.55 

0.00015 

33.079 

3.657 

4.697 

317.224 

1-7 

63-9-17-2-0 

42166.70 

0.00018 

33.043 

12.694 

0.581 

317.165 

1-8 

63-9-20-2-0 

42167.42 

0.00018 

33.010 

359.970 

16.282 

317.098 

1-9 

63-9-27-2-0 

42167.51 

0.00022 

33.046 

38.922 

344.162 

316.996 

1-10 

63-10-1-2-0 

42168.88 

0.00024 

33.039 

26.615 

0.433 

316.944 

1-11 

63-10-8-2-0 

42169.14 

0.00020 

33.013 

42.889 

350.866 

310.780 

1-12 

63-10-14-2-0 

42169.78 

0.00028 

32.982 

36.727 ; 

2.673 

316.813 

1-13 ! 

63-10-22-2-0 

42171.51 

0.00026 

32.993 

62.833 

344.246 

316.603 

1-14 ; 

63-10-30-0-0 

42171.09 

0.00028 

32.948 

29.865 

354.548 

316.570 

1-15 

63-11-6-0-0 

42172.15 

0.00025 

32.952 

36.699 

354.313 

316.328 

1-16 

63-11-12-5-0 | 

42172.51 i 

0.00031 ! 

32.920 

108.239 

3.425 

316.308 

2-1 

63-11-28-1-0 I 

42165.89 | 

0.00005 

32.920 

222.170 

203.901 

315.976 

2-2 

63-12-4-0-0 

42167.20 , 

0.00009 

32.892 

39.435 

17.564 

315.919 

2-3 

63-12-10-0-0 

42167.18 ! 

0.00010 1 

32.881 

51.942 

10.958 

315.877 

2-4 

63-12-16-17-0 

42168.17 

0.00007 | 

32.872 

300.000 

24.505 

315.735 

2-5 

64-1-6-17-0 

42168.01 

0.00013 

32.867 

332.997 

11.625 

315.544 

2-6 

64-1-9-6-0 

42169.90 

0.00015 

32.857 

165.031 

16.992 

315.469 

2-7 

64-1-20-21-0 

42171.43 

0.00012 

32.826 

29.098 

28.842 

315.300 

2-8 

64-1-29-20-0 

42171.91 

0.00019 

32.859 

37.956 

13.171 

315.212 

2-9 

64-2-5-16-0 

42172.89 

0.00019 

32.800 

321.168 

36.275 

315.075 

2-10 

64-2-10-19-0 

42173.31 

0.00014 

32.833 

32.517 

14.553 

314.982 

2-11 

64-2-17-17-0 

42174.89 

0.00019 

32.762 

347.774 

35.551 

314.883 
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Figure D1 — A portion of the celestial sphere at the epoch of orbit 2-5. 


On 6.0 January 1964, the hour angle of the vernal equinox west of Greenwich (expressed in hours, 
with 24 hours = 360 degrees) was 

6 hr 58 min 27.484 sec (from the Nautical Almanac) . 

On 7.0 January 1964, the hour angle of Q was 


7 hr 2 min 24.036 sec . 

Interpolating, the hour angle of Q on 6 January at 17 hours Universal Time was 

0 hr 1 min 15.042 sec, or 

0. 313 degrees west of Greenwich . 

In Figure Dl, the orbit angle 27.003 degrees is taken directly as 360 degrees minus the mean 
anomaly, because the orbit is nearly circular. The reported period for this orbit was 

T p = 1436.21696 min . 

The earth* s sidereal rotation period is taken to be 

T = 1436.06817 min . 

earth 

Thus, if the satellite is assumed to traverse orbit 2-5 at a nearly uniform rate, it will reach the 
celestial equator at a time when the Greenwich meridian has proceeded eastward from the epoch 

15.378 x 1436.21696/1436.06817 = 15.380 degrees . 
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Thus, the estimated geographic longitude of the ascending equator crossing nearest to the epoch 
of orbit 2-5 is 

Ascending equator crossing longitude = -(44.456 + 0.313 + 15.380) 

= - 60. 149 degrees (see Table 1). 

The estimated time of this crossing is 

15.380° 

15 °/}^ - 1.025 hours after the epoch . 

The crossing time (Table 1 in text) is thus estimated to be at 


6.751 January 1964 (18.025/24+6.0 January 1964) . 
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Appendix E 


Derivation of Exact Elliptic Integral of Drift Motion for 
a 24-Hour Near-Circular Orbit Satellite: Comparison 
of Exact Solution with Approximate Solutions 
for Periods Very Close to Synchronous 


The differential text equation 33 of 24-hour satellite drift is 
analogous to the equation describing the large-angle oscillations 
of a mathematical pendulum (see Reference 10), as in Figure E-l. 

The equation of angular motion of the mass m under the con- 
stant gravity force mg 0 is 

md 2 ( 16 ) _ - 

F t - mg 0 sin 9 - — ^ - m 10 ■ (El) 

Equation El can be rewritten as 

8 + (g o/l) sine = 0 . (E2) 



From the theory developed in Reference 10 (pp. 327-335), Equa- Fi^re El— Configuration of 

tion E2 has an integral 0 ' Whematical P endulom " 


t (time from 6 - 0) - (Z/g 0 ) 1/2 F(k, > (E3) 

where F(k, ±) is the elliptic integral of the first kind with argument (or amplitude) 



and modulus 

k = sin 0 mKl / 2 - 


Equation 33 in the text. 


y + A J2 sin 2y - 0 , 
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with maximum libration angle y Q9 can be put in the form of Equation E2 by the transformation of 
the dependent variable 


6 = 2 y 


(E4) 


with the parameter identification 


& 0 /l - 2A 


(E4a) 


Equation E4 implies the identification 


k - sin y Q , <£ - sin 1 (siny/siny 0 ) . 


(E4b) 


The pendulum solution (Equation E3), under the transformation E4 and identifications E4a and 
E4b, becomes 

t(time of drift libration from y = 0) = (l/2A 22 ) 1/2 Fjsiny 0 , ^sin" 1 siny/siny 0 )J . (E5) 

F(k , 0) , in its full integral form, is 

d<£ 

F = (1 -k 2 sin 2 ^) 1/2 ’ (E6 > 

Jo 


(where k 2 - sin 2 y 0 , sin 2 ± = sin 2 y/sin 2 y Q ) for the drift libration. In particular, when £ = tt/2, 
then y = y Q ; y 0 = 0; and the pendulum- drift libration has completed a quarter- period. 


Thus, from Equations E5 and E6, the full period of the long-term drift libration of the 24-hour 
satellite ground track about the nearest minor equatorial axis longitude is 

r 77 / 2 


T W 


(8/A 

Jo 


d jp 

(l - sin 2 y Q sin 2 ^) 1/2 


(E7) 


The adequacy of the Taylor series expansion approximation of the drift motion in the neighborhood 
of y 0 , given in Equation 50, may be tested against the exact drift solution implicit in Equation E5. 
Table El gives the evaluation of F for arguments within 5 degrees of > 0 = 60 degrees, using the 
integral tables in Reference 11. In Table El, Af is the change in the elliptic integral from the 
"stationary" configuration at y = 60 degrees or £ = 90 degrees; At = (l/2A 22 ) 1/2 AF; A 22 was com- 
puted from Equation 30a with the following gravity- earth constants and for the inclination of 
Syncom II: 


R 0 = 6378.2 km , 

a_ = 42166 km , 

S 
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J 22 = -1.7xlO-«, 

i - 33 degrees . 


7 ' gives the drift position as calculated from the first right-hand term of Equation 50 alone [the 
(At) 2 term]. y" gives the drift position as calculated from the first two right-hand terms of 50. 
The "actual" Syncom n drift in mid- August 1963 began, apparently, at a y 0 between 48 and 58 de- 
grees east of the minor axis. Thus, the 16 orbits chosen for the first drift period all should be 
well represented by the (At) 2 -only theory, within the rms error of the longitude observations. 
Similar exact calculations for y Q = 45, 50, and 55 degrees confirm the adequacy of the (At) 2 - 
only theory to apply to the second drift- period orbits. They also prove the contention in the sec- 
tion on approximations to exact drift solutions (p. 16 ff.) that, for reasonably small excursions 
from "synchronism", the convergence of the Taylor series (Equation 50) is adequate if additional 
terms are included only when they become of a certain minimum significance to the total drift 


Table El 


Exact and Approximate Drifts of a 24-Hour Satellite from a Stationary Configuration 
60° East of the Earth 1 s Minor Equatorial Axis 


y Q = 60 degrees 

7 

(degrees) 

7 } 

(degrees) 

y” 

(degrees) 

& 

(degrees) 

F 

(rad) 

(rad) 

At 

(days from) 

7 - 60° 

A 22 = 23.2 x 10” 6 rad/day 2 

60.0 

60.0 

60 

90 

2.1565 

- 

- 


59.0 

59.003 

59.000 

81.7967 

1.8730 

0.2835 

41.619 


58.0 

58.014 

58.001 

78.3056 

1.7564 

0.4001 

58.737 


57.0 

57.029 1 

56.999 

75.5595 

1.6671 

0.4894 

71.846 


56.0 

56.051 

56.999 

73.1938 | 

1.5923 

0.5642 

82.827 


55.0 

55.077 

54.996 

71.0617 

1.5265 

0.6300 

92.487 
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Appendix F 


List of Symbols 

A 22 Driving function causing drift and orbit expansion of a 24-hour satellite in a 
’’triaxial" earth-gravity field; a constant for a given 24-hour orbit inclination 

a, a s Instantaneous semimajor axis, and "momentarily synchronous" semimajor axis, of 
orbit of the 24-hour earth satellite. (a s , estimated to within 3 km, is 42166 km) 

a 0 , b 0 Major and minor equatorial radii of the "triaxial" earth [R 0 = (a 0 + b 0 )/2, accord- 
ing to definition in Wagner, op. cit.; see footnote, p. 7] 

a, (a - a s ) / /a s ; a nondimensional semimajor axis change for the 24- hour satellite's 
orbit, with respect to "momentarily synchronous" semimajor axis 

ds A small arc length of a space trajectory 
( d n ) • ( e n ) Determinable coefficients in drift and orbit- expansion equations 67 and 68 

F A gravity force per unit mass acting on a 24-hour satellite 
F( i ) Inclination factor in the triaxial driving function A 22 
F(k, £) Elliptic integral of the first kind with argument (or amplitude) ± and modulus k 

g 0 , g s Radial acceleration of the earth's gravity field at the earth’s surface and at alti- 
tude of the "synchronous" satellite 

i Inclination of orbit of the 24 -hour satellite 
J nm , X nm Spherical harmonic constants (order n, power m )* of the earth’s gravity potential 

M Mean anomaly of the satellite in its orbit: orbit angle (from center of the earth) 
from perigee to a point M in the orbit, where m = 27rt/T p , t being real time since 
perigee passage and T p period of the satellite's orbit 

m A test mass 

R 0 Mean equatorial radius of the earth (6378.2 km) 

*In the literature of spherical harmonics n is often called "degree,” and m "order." However, in speaking of the order of influence of 
the gravitational harmonics, n is generally used. For this reason, the author prefers the nomenclature "order” for n and "power” or 
"degree” for m when referring to gravitational harmonics. 
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T i T 

P s 


< >0 


( ) 


n 


C). ("). ( ) [n 3 
A( ) 

7 




e 

K r, <p 


K 


o 


Vk 




n 


CO 


e 


Orbital period for a satellite, and "momentarily synchronous" period of a 24- 
hour satellite (i.e., the earth T s sidereal rotation period) 

Time of "synchronism" from an arbitrary base time of reckoning T 
Real time 

Gravity potential of the earth 

Argument of perigee in a satellite orbit: orbit angle (from center of the earth) 
from ascending node to perigee 

Argument ( ) at start of dynamics under consideration 

Argument ( ) at a specified location n [except in Appendix C: ( )s; argument 

for the simulated trajectory] 

d ^ , d time differentials 

’ dt 2 dt n 
A small argument ( ) 

Geographic longitude (positive to the east) of the 24-hour satellite, or ascending 
equator crossing of the satellite's orbit with respect to longitude of the earth's 
minor equatorial axis 

Geographic longitude (positive to the east) of ascending equator crossing of the 
24-hour satellite's orbit with respect to the earth's minor equatorial axis' longi- 
tude location, at start of dynamics under consideration 

Argument from ascending node to satellite position for the 24-hour orbit 

Geographic longitude, geocentric radius, and geocentric latitude of the 24-hour 
satellite position 

"Initial" geographic longitude of the satellite, or ascending node of the 24-hour 
satellite's orbit at start of dynamics under consideration 

Geographic longitude difference between two "momentarily synchronous" 24- 
hour satellite configurations 

Earth’s gaussian gravity constant (3.986 x 10 5 km 3 /sec 2 ) 

Longitude location of the vernal equinox 

Earth’s sidereal rotation rate (0.7292115 X 10“ 4 rad/sec) 
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